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Abstract 

The quasinormal modes (QNM's) of gravitational systems modeled by 

the Klein-Gordon equation with effective potentials are studied in analogy 

to the QNM's of optical cavities. Conditions are given for the QNM's to 

form a complete set, i.e., for the Green's function to be expressible as a sum 

over QNM's, answering a conjecture by Price and Husain [Phys. Rev. Lett. 

68, 1973 (1992)]. In the cases where the QNM sum is divergent, procedures 

for regularization are given. The crucial condition for completeness is the 

existence of spatial discontinuities in the system, e.g., the discontinuity at the 

stellar surface in the model of Price and Husain. 
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Radiation from optical cavities is often analysed in terms of the "modes" of the cavity 
Because the waves escape, causing the total energy in the cavity to decrease, these 
"modes" are quasi-normal modes (QNM's). The observation of the QNM's from the outside 
immediately gives information on the spatial structure of the cavity (but not the spatial 
structure of the source of radiation inside the cavity), e.g., in obvious notation, uij ~ j-nc/L, 
where L is the length of a simple 1-dimensional optical cavity. 

Similarly, gravitational radiations from relativistic systems can be described by QNM's. 
Although "cavities" in gravitational systems, i.e., regions of spacetimes which reflect and 
scatter waves significantly, tend to be more leaky than optical cavities, we argued that QNM 
analysis and the optical analogy can still be very useful. Indeed the QNM's of black holes 
and relativistic stars have been subjects of much study. In numerical simulations, it is 
often found that the radiation observed in many black hole processes is dominated by 
the QNM's of the hole. For processes as violent and non-linear as the head-on collision of two 
black holes, the domination of the radiation by QNM's may seem surprising 0. However, it 
is easy to understand with the optical analogy — the distant observer sees only the QNM's 
of the optical cavity, but not the details of the radiation generation mechanism. 

Gravitational QNM radiation may be observed by LIGO and VIRGO |^ in the next 
decade, and may therefore reveal spacetime structures of various gravitational systems, e.g., 
the strong field region around a black hole — in much the same way as the spectrum of a 
laser permits a distant observer to infer some characteristics of the cavity. 

This exciting possibility calls for the study of the general properties of QNM's of gravi- 
tational systems. For the interpretation and extraction of observational data, one would like 
to know, for example, how the QNM frequencies of a black hole are perturbed by a massive 
accretion disk around it. 

A fundamental question, raised by Price and Husain is whether QNM's can form 
complete sets. This question is not only of theoretical interest, but also of technical impor- 
tance, e.g., the application of the usual Rayleigh perturbation theory to calculate frequency 
shifts hinges on the existence of a complete set. The answer to this question is thought to 
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be in general negative, as the system giving rise to QNM's is nonhermitian. Yet Price and 
Husain ^ have given a model of relativistic stellar oscillations that does have a complete 
set of QNM's. 

In this Letter, we consider this question in a general context, and give conditions for 
completeness to hold. The crucial condition is the existence of spatial discontinuities. (In 
the model of 0, this discontinuity corresponds to the stellar surface Such a condition 
is not surprising: a discrete sum of QNM's can at best be complete over a finite interval; the 
discontinuities provide a natural demarcation of this interval, analogous to the boundaries 
of an optical cavity. 

The propagation of waves in curved space is often modeled by the Klein-Gordon (KG) 
equation 

D<Pix, t) ^ [d^ -dl + V{x)\ 0(x, t) = , (1) 

with the outgoing wave boundary condition at infinity. The potential V{x), assumed to 
be positive, bounded and vanishing at infinity, describes the scattering of wave by the 
background geometry. We show that the QNM's of (1) form a complete set, in the sense 
that the Green's function can be expressed as a sum over QNM's (cf. (10) below), provided 
the following three conditions hold: (i) V{x) is everywhere finite, and vanishes sufficiently 
rapidly, as x ^ cxo, in a sense to be defined below; (ii) there are spatial discontinuities 
demarcating a finite interval; (iii) consideration is limited to certain domains of spacetime. 
In particular, for propagation from a source point y to an observation point x in a time t, 
the QNM's give a complete description for (a) y inside the interval, x outside the interval, 
and if: > tp(x, y), for a certain tp\ or (b) both y and x inside the interval, and t > (with the 
retarded Green's function being zero for t < ). Case (b) is analogous to the completeness 
of normal modes in hermitian systems defined on a finite interval. 

The above conditions allow dispersion, backscatter, as well as differences in the damping 
times, the absence of which has been conjectured to be important []^. The results are valid 
to all orders in the rate of leakage. (Lowest order results would be trivial, since the system 



becomes hermitian in that limit.) These results extend our previous work on the wave 
equation |I0| . 

The spatial coordinate x in (1) often represents a radial variable, so first consider a 
half-line problem (x > 0) with (f){x = 0,t) = 0. The QNM's are eigen-solutions to the 
time-independent KG equation 

_ Q2 ^ = , (2) 

with the boundary conditions 0(0) = and the outgoing wave condition at infinity. The 
eigenfunctions and eigenvalues are denoted as 0(x) = fj{x) and uj = ojj. 

The retarded Green's function for the system is defined by DG{x, y; t) = 6{x — y)S{t) 
with G = for t < 0, and the boundary conditions are (i) G = for either x = or 
y = 0, and (ii) the outgoing wave condition as either x —>■ oo or y oo. In the frequency 
domain, and henceforth choosing y < x, G{x,y;Lj) = f{uj,y)g{uj,x)/W{u), where / and g 
are solutions to (2) with the boundary conditions /(cj, x = 0) = 0; /'(cu, x = 0) = 1 [|11[], and 
g{uj,x) —>■ exp{iujx) as x ^ oo, where ' = d/dx. 

The Wronskian W{uj) = gf — fg' is independent of x, and its zeros are the QNM 
frequencies ojj, with f{ujj,x) oc g{Ljj,x) oc fj{x). For simple zeros |T2[, the residues of 



G{x,y]Lj) are Kj = f{Ljj,y)g{ujj,x)/[dW{uj = ujj)/duj]. Multiple zeros can be handled 
readily. To analyse the denominator dW/du, start with the defining equation for f{u!j,x) 
and g{u,x). The usual manipulations lead to 

rX 

dx fijjJj, x)g{u!, x) 

= [gi^, x)f'{ujj, x) - g'{uj, x)f{ujj, x)] J /{uj^ - uj]) , (3) 

where the integral is taken along any contour from x = to x = X. Use the outgoing wave 
condition to evaluate / and g at the upper limit, and take uj —>■ ujj by I'Hospital's rule. This 
then gives 

rX 

/ dx f{uj,x)g{ujj,x) 
Jo 

= _ [51^/5^ + ,/^] /(2^^.) , (4) 



for X oo. Since / and g are proportional at these poles, a generalized norm of the QNM's 
can be defined as 

« I » = lim r dx f,{xf + if,{Xf/ {2u;,) , (5) 

A— >oo Jo 

with which the residue of G at ujj can be expressed as 

= -mm/ [2^. « /. I /. »] • (6) 
The generalized norm, introduced in other contexts |]l],|l3|,|l4| , involves rather than 



and is therefore complex. Although each term on the right of (5) does not have a limit as 
X oo, the limit exists for the sum. 

Next write G in terms of G by a Fourier integral, and distort the contour to a large 
semicircle in the lower half u plane. One then sees that 

Z ■ UJj ^ Jj I Jj » 

The sum comes from the zeros of W{u!), Ic is the integral along a semicircle at infinity, and 
Is comes from any singularities of / and g in the lower half u plane. It remains to show 
that (i) / and g are regular in u, in which case Is = 0; and (ii) Ic vanishes if there is a 
discontinuity in V{x) at some x = a > 0. 

The first step is straightforward, by appealing to well-known results in the quantum 



theory of scattering [|T^]. Since / and g satisfy the Schroedinger equation with energy , 



they are analytic functions of u if the potential is bounded and "has no tail" [0, in the 
sense that dx x \ V{x) \ < oo, and dx x e"^ | V{x) \ < oo for any a > . Our 
work on potentials that have a tail, e.g., inverse power laws, will be reported elsewhere. 

Secondly, on the large semicircle G can be estimated by the WKB approximation. 
As |u;| oo WKB approximation is valid except at points where the potential has dis- 
continuities, which can be handled by connecting the WKB solutions across them. Let 
0(x) = exp[iS{x)] be a solution of (2); then S{x) ~ ±/ k{x)dx, where k(x) = [cu'^ — V{x)Y^'^. 
G can then be obtained in terms of (f){x) and k{x). 



Consider a potential with a step discontinuity at x = a. The reflection coefficient is 
R = [S'{a^) — 5"(a+)]/[S"(a^)* + S"(a+)] ~ uj~^ at high frequencies. For simphcity we give 
the argument only for Case b: < y < x < a. It is straightforward to show that 

~ _ sin[/(0,?/)][e-*-^(^'") + i?e^^(^''^)] 
^'^'^ ~ ^A:(a;)A:(y)[e-*^(M + i?e^^(o.«)] 

where I{u, v) = /J k{x)dx ~ uj{v — u). Now on the semicircle uj = ujr + iuJi = Ce*^, it < 6 < 
27T, as uj = CsinO — oo, both the numerator and the denominator of G are dominated 
by the term proportional to R, and 

As C — > oo, this vanishes for t > 0. This conclusion remains valid if V{x) has a discontinuity 
only in its p-th derivative {p = 0, 1,2, ...), since the reflection coefficient R would vary as 
Lj-(p+2) ng Thus we have proved (apart from a technical hitch mentioned below) that for a 



discontinuous potential, the QNM's are complete for t > 0. For Case a, i.e. 0<y<a<a;, a 
similar estimate gives, in place of (9), e~*'^^*~^~^"'"^"V(i?c<j). This vanishes and consequently 
completeness holds only if t > tp{x, y) = max(x + y — 2a, 0). 

For problems on a full line — oo < x < — oo (as in the case of Schwarzschild hole), / 
is still defined as the solution satisfying the left boundary condition, which is now a unit 
outgoing wave as x ^ — oo; all the arguments remain the same. Generalization to multiple 
discontinuities is likewise straightforward. In this case, completeness in Case b holds in the 
interval between the leftmost and the rightmost discontinuities. 

However, there is one technical hitch. The QNM sum is in fact a divergent series, but 



converges to the correct answer when regularised in a standard way [|T7|- The simplest 
regularization is to ffist invoke symmetry and only keep modes with Re Uj > (hereafter 
denoted as j > 0) in the sum in (7); secondly replace all times t by t — ir, with t O'^ . This 
imphes the use of a regulating factor /j (r) = exp(— cUj-r) multiplying the contribution of the 
j-th QNM. The need for regularization can be seen as follows. The proof that the integrand 
vanishes on the large semicircle as C — > oo is not valid near the real axis, in particular, 
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where Im uj = O(logC). In this part, a factor such as exp(Im ucr) {a real and positive, e.g. 
a = t — X — y + 2a) behaves hke a power of C rather than exponentially, and it is necessary 
to examine the integrand more carefully ||I8[. If i? ~ cj^'', the powers of u) in the integrand 



for G are inadequate to control the contribution along the semicircle if g > 1 [0]. For the 
KG equation with V{x) discontinuous in its p-th derivative, q = p + 2, so regularization is 
always necessary. For the wave equation with a discontinuity in the p-th derivative of the 
refractive index n{x), q = p, so regularization is not necessary for p = 0. The difference by 
two powers can be understood with the transformation relating the wave equation to the 



KG equation ||T0 |. 

For large p, the divergence is more severe, the regulated sum converges more slowly to 
the correct result as r ^ 0, and the use of a small finite r becomes increasingly inaccurate. 
Thus, in practice, the QNM sum is only useful for relatively strong discontinuities, e.g. 
p = 0,l. 

To summarize, we have proved, under the conditions stated, that the Green's function 
G is expressible, in a certain domain of x, y, t, as the regulated QNM sum 



G = lim Re 



j>0 



(10) 



Considering G as t — > 0+, 



6{x-y)=G{x,y;t = 0+) 



lim Re 



(11) 



valid for x, y & [0, a]. This is perhaps the more familiar notion of completeness |]T9 



We now focus on (11), noting that such a resolution of the identity will be useful for a 
variety of problems similar to those of hermitian systems. The QNM representation of G 
for t > will be further discussed elsewhere in the context of initial value problems. As an 
example, consider a potential defined on the full line, with V{x) = Vi > for < x < a, and 
zero otherwise. This potential has discontinuities and no tail, thus satisfying the conditions 
necessary for our results. The QNM's of this system are easily found. We consider (11), 



but as an approximate equality for some small nonzero regularization parameter r, and in 
a distribution sense, e.g. 



ry2 

Re E fji^) / dyf,{y)I,{T)/ < /, | /, > 



j>0 

e{x-y^)-e{x-y2) , (12) 



for X, yi, y2 & [0, a]. Figure 1 shows that for small r, the partial sum indeed converges close 
to the right hand side. Without the regulator, the sum does not converge. 

This example further illustrates why the series needs to be regulated. Suppose that for 
a full-line problem, the reflection coefficients on both sides behave as R{u!) ~ uj^'^ at high 
frequencies {q = 2 for the present situation). By setting the denominator in (8) to zero ||20|| , 



it is straightforward to show that Im uja ~ — glog In the region < x,y < a, a typical 
term in the product fj{x)fj{y) goes as 

where A = max(|x + y — a\,\x ~ y\). Since the maximum value of A is a, without the 
regulator the worst behavior of the summand in (10) is The sum would not converge 

without a regulating factor if g > 1. The culprit is Im Uja ~ —q\og \j\, from which one 
can say unequivocally that the need for regulating the sum is an intrinsic property of KG 
open systems. 

We have shown that for potentials that have no tail but which contain spatial disconti- 
nuities, the QNM sum is complete, if suitably regulated. This result settles a question in 
the literature 0. The generic need for regularization, in particular, distinguishes the KG 
equation from the wave equation. When there is a tail, there will be scattering of the wave 
from asymptotic regions of space, leading to power-law type long time behavior, which will 
be discussed elsewhere. 

In so far as the Green's function G provides the solution to all the dynamics, the QNM 
expansion of G will lead to a variety of physical applications, much as the normal mode 
expansion of hermitian systems. 
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FIGURES 

FIG. 1. Smooth envelope forming upper bound of logi^E vs J/10^, where E is the difference 
between the right hand side of (12) and the partial sum on the left hand side, up to J terms. The 
parameters are Vi = 10^, a = 1, yi = 0, y2 = 0.5 and x = 0.25, r = 5 x 10~^ (curve 1); x = 0.25, 
r = 5 X 10-^ (curve 2); x = 0.75, r = 5 x IQ-^ (curve 3); x = 0.75, r = 5 x 10"^ (curve 4). 
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